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COMPENDIUM  OP  IMPEDANCE  FORMULAS 

ABSTRACT 

A  list  of  impedance  formulas  is  presented  which  can  be  applied  to  the 
vibration  excitation  of  simple  structures  such  as  beams/  plates,  rings, 
and  beam-plate  systems.  Infinite,  semi-infinite,  and  finite  systems 
are  considered.  Most  of  the  formulas  have  been  taken  from  the  exist¬ 
ing  literature;  others  are  derived.  In  the  latter  case  short  deri¬ 
vations  are  given  in  the  Appendixes. 
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LIST  OP  SYMBOLS 


bending  waves 
longitudinal  waves 
torsional  waves 

=  frequency 
=  thickness 

=  —  =  wave  number  for  bending  waves 
=  length 

=  surface  mass  or  linear  mass 
=  reflection  coefficient 
=  particle  velocity 
=  angular  velocity 
=  coordinates 

=  bending  stiffness  of  beams;  for  rectangular  cross  section 

r,  Eh3b 
B  -  ~12~ 

=  torsional  stiffness  of  beams 

=  bending  rigidity  of  plates;  for  rectangular  cross  section 


=  radius 

=  width  of  a  beam 
=  propagation  velocity  for 
=  propagation  velocity  for 
=  propagation  velocity  for 

=  transmission  coefficient 


D  = 


Eh3 

12(l-h2) 
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D' 

E 

F 

G 

L 

M 

P 


A 

H 

P 

co  =  2irf 


=  torsional  rigidity  of  plates 
=  Young's  modulus 
=  force 

=  shear  stiffness 
=  length 
=  moment 

=  mechanical  power 

=  reflection  coefficient  for  bending  wave  near  fields 
=  cross  sectional  area  of  beams 
=  force  impedance 
=  moment  impedance 
=  source  impedance 
=  loss  factor 
=  radius  of  gyration 
=  wavelength 
=  Poisson's  number 
=  density 

=  angular  frequency 
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I.  INTRODUCTION 

The  concept  of  mechanical  impedance  is  very  useful  in  solving  vi¬ 
bration  problems  because  it  allows  expression  of  the  energy  transfer 
from  a  vibration  source  to  a  structure,  and  from  one  structure  to 
another,  in  fairly  simple  terms.  It  is  an  especially  convenient  _____ 
concept  for  engineers  with  background  in  electrical  engineering  be¬ 
cause  of  the  analogy  with  electrical  and  mechanical  impedances. 

A  particular  example  of  the  usefulness  of  the  impedance  concept  is 

the  design  of  vibration  mounts,  where  the  vibration  reduction  can 

be  predicted  rather  accurately  if  all  the  impedances  involved  are 

known.  But  there  is  one  difficulty:  In  many  practical  cases  the 

impedances  are  not  known,  and  therefore  one  is  forced  to  represent 

certain  parts  of  structures  by  lumped  masses,  springs  and  dashpots. 

This  approach  is  convenient  and  useful  for  low  frequencies  or  for 

heavy  and  very  compact  structures,  but  for  higher  frequencies  and 

1-4/ 

light  structures  it  may  be  very  misleading. - '  Thus,  one  is  faced 

with  the  problem  of  getting  more  information  about  the  impedances 
of  real  structures. 

It  is  hoped  that  the  present  report  will  provide  some  of  this  infor¬ 
mation.  The  report  gives  the  impedance'  of  structures  which  are  more 
complic^ijed  than  lumped  masses  and  springs,  but  still- simple  enough 

to  be' treated  mathematically. 

» 

Most  of  the  formulas  of  this  report  deal  with  infinite  systems. 

This  is  not  as  serious  a  restriction  as  it  might  appear  provided 
that  one  is  not  interested  in  the  response  of  a  system  at  a  given 
frequency,  but  rather  in  the  average  behavior  in  a  frequency  band. 


/ 
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It  will  even  be  shown  in  Section  II-A4  that  the  input  impedance  of 
an  infinite  system  governs  the  flow  of  mechanical  power  into  a 
finite  system  of  the  same  kind. 


There  are  two  types  of  impedances  listed  in  this  report;  one  of  them 
is  the  force  impedance  Zp  which  is  given  by 


(1-1) 


(P  =  exciting  force,  v  =  particle  velocity  at  the  excitation  point); 
the  other  one  is  the  moment  impedance  given  by 


(1-2) 


(M  =  exciting  moment,  w  =  angular  velocity  at  the  excitation  point ) . 


In  both  cases  only  point  forces  or  moments  are  considered.  This 
means  that  the  excitation  is  localized  in  a  region  which  is  very 
small  compared  to  the  wavelength  in  the  structure,  otherwise  the 
formulas  cannot  be  applied. 


The  more  familiar  force  impedance  is  used  if  all  the  excitation  is 
given  by  a  force;  one  example  might  be  the  excitation  of  a  structure 
by  an  engine  (pump,  etc.)  provided  that  the  vibrations  of  the 
engine  (etc.)  are  purely  perpendicular  to  the  structure.  If  there 
is  also  some  "rocking"  of  the  engine,  the  exciting  moment,  and  there¬ 
fore  the  moment  impedance^,  has  to  be  taken  into  account,  too. 


5  * 
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Since  most  of  the  formulas  given  in  this  report  have  restricted 
ranges  of  applicability.,  an  -attempt  has  been  made  to  give  not  only 
the  equations  but  also  the  corresponding  restrictions.  With  regard 
to  sign  convention  we  emphasize  that  we  have  assumed  harmonic  motion 
of  angular  frequency  o>  with  the  time  dependence  expressed  by  e^Q^. 
This  is  in  agreement  with  Cremer's  work^/  from  which  many  of  the 
formulas  were  taken.  If  harmonic  motion  of  the  form  e-iayt  is 
assumed,  -i  instead  of  +i  must  be  inserted  in  the  formulas.'  Thus, 
in  the  present  report  an  impedance  with  a  positive  imaginary  part 
.is  mass-like,  while  an  impedance  with  a  negative  imaginary  part  is 
stiffness-like . 
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II.  FORCE  IMPEDANCE  OF  BEAMS 


Excitation  of  Bending  Waves 


The  impedance  of  beams  or  systems  of  beams  (for  sketches  of  the 
cases  considered  in  the  following,  see  Fig.  l)  which  are  excited  by 
a  point  force  acting  perpendicularly  to  the  beam  can  be  calculated 


by  solving  the  one -dimensional  bending  wave  equation 


“  4^v  =  0  (Il-l) 

dxH 


under  the  proper  boundary  conditions.  In  Eq.  (il-l),  v  is  the 
transverse  velocity  of  the  beam 


(II-2) 


the  bending  wave  number,  m  the  mass  per  length,  S  the  cross-sectional 
area  of  the  beam,  k  the  radius  of  gyration,  and  E  is  Young's  modulus. 

Equation  (il-l)  is  valid,  provided  that  the  bending  wavelength  A  = 
27r/k  Is  much  larger  than  the  thickness  h  of  the  beam.  In  general, 

Eq.  (ZI-l)  can  be  used  if 

A  >  6h  .  (II-3) 

A  similar,  somewhat  less  stringent  restriction  holds  for  the  width 
of  the  beam. 
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1 .  Impedance  of  Uniform  Infinite  Beams 

The  term  "infinite  beam"  does  not  mean  that  the  beam  must  actually 
be  of  infinite  length.  All  the  following  results  are  correct  if  the 
beam  is  highly  damped  or  otherwise-  terminated  in  such  a  manner  so 
that  no  reflected  bending  waves  come  back  from  that  side  of  the  beam 
which  is  assumed  to  go'  to  infinity.  Methods  for  terminating  beams 
have  been  described  by  Kurtze,  Tamm,  and  Vogel  & 

a)  Beam  extending  from  x  =  -»  to  x  =  -H° 

In  this  case  (investigated  by  Cremer^  the  boundary  conditions  at 
the  driving  point  are,  first  that  the  angular  velocity  is  zero,  and  • 
second, .that  the  exciting  force  is  equal  to  the  transverse  force  of 
the  bending  wave.  This  gives  for  the  impedance 


=  2mc(l+i) 


2SEfC?k3 


(1+i) 


(II-4) 


where  c  =  =  bending  wave  velocity.  The  geometry  for  this  case  is 

shown  in  Fig.  la. 

b)  Beam  extending  from  x  =  0  to  ,x  =  °° 

This  case  was  also  considered  by  Cremer.-^  See  Fig.  lb  for  a  sketch 
of  the  geometry.  Using  as  boundary  conditions,  first  that  the  excit¬ 
ing  force  is  equal  to  the  transverse  force  of  the  bending  wave  and 
second,  that  there  is  no  bending  moment,  he  obtained 


ZF  " 


(l+i ) 


(H-5) 
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2.  Impedance  of  Beams  with  One  Reflecting  Device 

To  make  the  following  formulas  as  general  as  possible  we  define  a 
reflecting  device  as  any  inhomogeneity  in  material  or  cross-section 
which  causes  a  reflection  of  bending  waves .  A  list  of  reflection 
coefficients  r  for  some  special  cases  is  given  in  Appendix  A.  This 
list  includes  the  reflection  coefficients  for  the  reflected  non¬ 
propagating  wave  (near  field)  as  well  as  the  reflected  propagating 
wave  (far  field).  We  shal-1  always  assume  that  the  distance  between 
the  driving  force  and  the  closest  reflecting  device  is  greater  than 
half  a  wavelength.  This  allows  us  to  neglect  the  complications 
caused  by  the  reflected  non-propagating  wave  (near  field). 

.a)  Beams  extending  from  x  =  to  x  =  -p  with  a  reflecting 
device  at  x  =  -1;  (1  may  be  positive  or  negative) 

The  impedance  of  this  structure  (shown  in  Pig.  lc)  can  be  found  as  a 
limiting  case  of  Eq.  (B6)  in  Appendix  B.  This  equation  gives  for 
vz  =  rv  rL  =  0 


4mc 


1-i  +  re 


-2ikii 


(II-6) 


b)  Beams  extending  from  x  =  0  to  x  =  -F°  with  a  reflecting 
device  at  x  =  l 

In  this  case  (shown  in  Pig.  Id)  a  derivation  very  similar  to  the  one 
in  Appendix  B  can  be  made .  The  main  difference  is  that  the  impedance 
of  the  initial  wave  is  given  by  Eq.  (II-5)  instead  of  Eq.  (II-4) . 
Furthermore,  use  is  made  of  the  fact  that  for  a  free  end  (at  x  =  0) 
the  reflection  coefficients  are  -i  and  1-i,  respectively.  Thus,  we 
get  the  following  expression  for  the  impedance 
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(1+1) 


1  +  lre-21ki 
1  +  riT2i“ 


■  (11-7) 


For  many  applications  it  might  be  of  interest  to  find  the  mean  square 

- p- 

velocity  vQ  at  the  driving  point  if  the  beam  is  excited  by  a  broad 
frequency  band.  Contour  integration  of  Eq.  (II-7)  shows  that- 


2  |f1  2  1  +  I  r  |  2 
m2c2  1  -  r|2 


(II “8) 


This  result  can  be  extended  easily  to -beams  with  internal  losses 
characterized  by  a  loss  factor  q.  .In  this  case  we  have 


2  1  F 1 


2  2 
c  m 


1 _ + 

*1  - 


(11-9) 


3 .  Impedance  of  Beams  with  More  Reflecting  Devices 

a)  Beams  extending  from  x  =  -°°  to  x  =  4-°°,  with  reflecting 
devices  at  x  -  -L  and  x  =  l  (L  and  l  positive) 

This  case,  shown  in  Fig.  le,  is  treated  in  the  second  part  of. 
Appendix  B.  Eq.  (B6)  derived  therein  leads  to  the  following  ex¬ 
pression  for  the  force  impedance 


-  4mc 
F 


1-ririe- 


-Sik(l+L) 


1_i+rLe-2ikL+^e-2ik.+(l+i)r^Le-2ik(i;+L) 


(11-10) 


(For  rg  =0,  Eq.  (II-10)  is  identical  with  Eq.  (II-6)). 
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The  mean  square  velocity  at  the  driving  point  for  broadband 
exci-tation  is  obtained  by  integrating  Eq.  (b6)  over  frequency. 
This  gives 


2  2  + 


1(-2  2 
16c  m 


aa_Li? 

l'  12 


1  “  ,rLri' 


(H-11) 


If  the  damping  of  the  beam  is  significant,  Eq.  (il-ll)  can  be  ex¬ 
tended  by  inserting  |r^rL|2  e-Rk(i+L)  in  place  of  |rgrL|2. 

b )  Beams  extending  from  x  =  -”  to  x  =  +°°,  with  several 
reflecting  devices  on  both  sides  of  the  driving  point. 

Equation  (11-10)  also  holds  for  this  case.  We  only  have  to  assume 
that  instead  of-  simple  reflecting  devices  at  x  =  -L  and  x' =  H, 
there  are  t-wo  "black  boxes"  containing  more  than  one  inhomogeneity. 
The  problem  is  now  reduced  to  determining  the  reflection  coeffi¬ 
cients  r^  and  r2  of  the  "black  boxes."  Unfortunately  the  reflection 
coefficients  ’of  systems  more  complicated  than  those  given  in 
Appendix  A  are  hard  to  compute.  But  the  average  behavior,  which  is 
given  by  Eq.  (il-ll),  can  easily  be  found  if  we  remember  that  (if 
no  internal  losses  occur)  the  energy  coming  to  a  reflecting  point 
is  equal  to  the  sum  of  the  reflected  and  transmitted  energy.  Or 
in  terms  of  the  reflection  coefficient  |r|  and  the  transmission 
coefficient  d 

| r | 2  =  1  -  | d | 2  .  (11-12) 

If  we  insert  this  equation  in  (il-ll),  we  get 
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“fjf.al 


I  2 


l6c2m2 


3  -  -  |d2| 

p  p  2 

|d2Hd2|-|d1d2| 


l6c2m2  |d2|+|d2 


(II-13) 


The  latter  approximation  holds  for  d^  «  1,  d2  «  1 


Equation  (±1-13)  is  correct  if  the  damping  is  small.  This  condition 
can  be  written  as: 


j d2 1  +  | d2 1 2  >  e~rik('g+L)  . 

4.  The  Power  Flow  Into  Finite  Beams 

At  this  stage  it  seems  appropriate  to  give  a  general  expression  for 
the  average  mechanical  power  that  flows  into  a  beam  from  a  vibration 
source.  To  this  end  we  consider  a  finite  beam  of  length  Z  and  loss 
factor  tj o  This  beam  is  excited  by  a  point  force  F  at  the  point  ZQ. 
The  velocity  of  the  beam  can  then  be  found  in  terms  of  the  eigen¬ 
functions  #n(x)  and  eigenfrequencles  can.  Thus  we  get  (for  example 
see  Ref.  7} 


oo 

.z  Y 


iw»n(i0) 


m  nio  m2  -m2  +  iT)m2 


Therefore  the  mechanical  power  P 


P  =  |  Re 


>  = 


2m 

n=0 


ax 

transmitted  into  the  beam  is 

^n  ^n  (foj 

go2  -  cd2  +  iiia)2|2  /02 


(11-14) 


(II-15) 
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Averaging  over  all  possible  exciting  points  £  and  over  a  broad 
frequency  band  Ao  gives 


p'  _  7T  F2  AN 
7  mJ  Acd 


(11-16) 


Here  AN  is  the  number  of  resonance  frequencies  in  the  frequency 
band  Aon  If  there  are  several  resonance  frequencies  within  the 
band  of  interest,  we  can  approximate  Eq.  (II-16)  by  inserting  the 
asymptotic  formula 


AN  ~  c£_ 
Aco  ~  2tt 


(11-17) 


Thus  we  get 


P 


a  F 


2  1 
8cm 


(11-18) 


.which  is  exactly  the  formula  that  gives  the  power  that  Is  trans¬ 
mitted  from  a  point  source  into  an  infinite  beam  (see  Eq.  (II-4)). 


Equations  (II-16)  and  (II-18)  show  that  the  average  power  is  in¬ 
dependent  of  the  loss  factor  q  and  the  length  of  the  beam.  This 
fact  has  an  interesting  consequence  for  the  power  transmitted  into 
a  beam  with  reflecting  devices. 


If  the  reflection  coefficients  are  rather  high  it  is  a  good  approxi¬ 
mation  to  assume  that  the  part  of  the  beam  which  is  driven  and  which 
is  between  two  reflecting  devices  (which  are  a  distance  £  apart) 
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behaves  like  a  finite  beam  of  length  i  with  additional  damping. 
This  means  that  Eqs .  (II-16)  and  (II-18)  hold  in  this  case,  too; 
or  to  put  it- in  another  .way,  the  average  mechanical  power  which 
is  transmitted  into  a  beam  is,  approximately,  unaffected  by  the 
presence  of  reflecting  devices,  and  can  easily  be  computed  if 
the  impedance  of  the  corresponding  infinite  system  is  known. 


5 .  Impedance. of  Beams  Driven  with  Non-Zero  Impedance  Sources 

In  all  formulas  given  above  we  have  assumed  that  the  source  does 
not  influence  the  beam  impedance;  i.e.,  the  source  impedance  Z0 

S  " 

has  been  assumed  zero,  or  at  least  much  smaller  than  the  impedance 

of  the  beam.  But  there  are  many  cases  in  practice  where  the  source 

impedance  Z  cannot  be  neglected  (e.g.,  a  heavy  machine  mounted  on 
s 

a  beam) . 


For  the  cases  Gee. (Il-Al)  (without  reflecting  devices)  the  in¬ 
fluence  of  the  source  can  easily  be  computed  by  adding  the  source 
impedance  to  the  beam  impedance  as  given  in  Eqs.  (II-4)  or  (II-5). 


For  all  other  cases  (with  reflecting  devices)  simply  adding  the 

source  impedance  would  be  incorrect,  because  of  the  interaction  at 

the  driving  point  of  the  source  and  the  waves  reflected  back  from 

the  reflecting  devices.  We  therefore  must  know  the  reflection 

coefficients  for  the  propagating  waves  and  non-propagating  waves 

at  the  driving  point.  If  these  quantities  are  denoted  r  and  R 

s  s 

respectively,  we  get  by  extension  of  Section  II. A. 2. a 


2mc (l+i) 


(l-l)(l-rr,  '-***) 
l-i+re~2ilce  (1+R  +ir  ) 

5  0 


(11-20 
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Extension  of  Eq\  (II-18)  in  Section  II.A.2.b  gives  [ 


Z„  = 


+  z 

2  + 


r  re 
s 


-2iki 


1+1  re"2iki(l+Rs-rs) 


(11-21) 


We  see  that  Eqs.  (11-20)  and  (11-21 )  cannot  he  obtained  by  adding 
Zg  to  Eq.  (II-6)  or  Eq.  (II-7).  The  reason  for  this  somewhat  sur¬ 
prising  result  is  that  in  general  the  reflection  coefficient  r  at 

s 

the  source  does  not  depend  only  on  the  force  impedance  Z  of  the 

s 


source  but  also  on  the  moment  impedance  at  the  driving  point, 
can  be  shown  for  example  that  Eq.  (II-21)  becomes 


Z  =  ~(l+i) 


1+ire 


-2iki 


1+re 


-2iki  +  Zs 


It 


(11-22) 


(compare  with  Eq.  II-7)  provided  that  the  moment  impedance  of  the 
source  vanishes.  Equation  (11-22)  can  be  proved  by  putting  v  =  0 
in  Eqs.  (A13)  and  (Al4)  and  introducing  the  resulting  expressions 
into  Eq.  (11-21 ).  . 

B.  Excitation  of  Longitudinal  Waves 

The  impedance  of  beams  (for  sketches  of  the  geometries  see  Fig.  2) 
excited  by  a  force  acting  in  the  direction  of  the  beam  axis  can  be 
computed  from  the  equation  for  longitudinal  vibrations  in  beams 


d2v 

dxS 


v  =  0 


(11-23) 


In  this  equation  v  is  the  longitudinal  particle  velocity. 
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(11-24) 


is  the  velocity  of  propagation  for  longitudinal  waves,  E  is 'Young's 
modulus,  and  p  is  the  density  of  the  beam  material. 

Equations  (11-23)  and  (11-24)  are  valid  only  when  the  cross-sectional 
dimensions  of  the  beam  "are  very  small  compared  to  the  longitudinal 
wavelength. 

1 .  Impedance  of  Uniform  Infinite  Beams 

The  longitudinal  impedance  of  a  beam  extending  from  x  =  -»  to* 
x  =  +oo  and  driven  at  some  point  in  between  is  given  by  (see  Fig.  2a) 

Zp  =  2S  pcL  =  2SyEp  (II-25) 

(S  =  cross  section  of  the  beam). 

x  =  +o>  and  driven  at  x  =  0  the 

(11-26) 

2.  Impedance  of  Beams  with  Reflecting  Devices 

The  impedance  of  beams  with  discontinuities  can  be  calculated  by 
using  the  methods  given  in  Appendix  B.  Since  all  waves  are  propa¬ 
gating,  the  formulas  obtained  are  correct  even  if  the  distance  be¬ 
tween  the  source  and  the  discontinuity  is  much  smaller  than  a  wave¬ 
length.  Furthermore  there  is  a  simple  relation  between  the  reflection 
coefficient  and  &he  Impedance  looking  into  the  discontinuity;  thus  the 
behavior  of  the  beam  can  be  expressed  solely  by  impedances  in  a  fairly 
simple  way. 


For  a  beam  extending  from  x  =  0  to 
impedance  is  (see  Fig.  2b) 
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a)  Beams  extending  from  x  =  0  to  x  =  °°  with  a  reflecting 
device  at  x  =  1  (i  >  0) 

If  the  reflection  coefficient  at  x  =  1  Is  r  we  get  (see  Pig.  2c) 


1-re 


-2i\a 


1+re 


-2ik,2 


(H-27) 


If  Zjj  is  the  impedance  looking  into  the  discontinuity  the  reflection 
coefficient  can  be  written  as 


ZD  -  SVEP 

zD  +  syE?  • 


(11-28) 


b)  Beams  extending  from  x  =  -»  to  x  =  with  reflecting 
■  devices  at  x  =  I  and  x  =  -L 


In  this  case  the  impedance  turns  out  to  be  (see  Fig.  2d) 

,-2ikU+L) 

>-2ik(i+L) 


Zp  =  2Sy^ 


1  "  rirL  e 


1  +  r£e 


-2iki  ;  -2ikL  +-r^Te 


(II-29) 


(r^  =  reflection  coefficient  at  x  =  &;  r^  =  reflection  coefficient  at 
x  =  -L )  , 

The  values  for  r^  and  r^  can  be  obtained  from  Eq.  ( 11-28)  if  the 
impedances  at  x  =  Z  and  x  =  -L  are  known. 
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Finally,  it  should  be  mentioned  that  for  longitudinal  waves  the 
effect  of  a  finite  source  impedance  Zg  can  always  be  taken  into  • 
account  simply  by  adding  it  to  the  impedance  of  the  structure.* 
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III.  FORCE  IMPEDANCE  OF  PLATES 
A.  Impedance  of  Infinite  Isotropic  Plates 

If  a  force  F  acts  perpendicularly  to  a  plate,  the  impedance  can 
he  found  by  solving  the  plate  equation  which  is  an  extension  of 
Eq.  (ll-l).  If  A  is  the  Laplacian  operator,  the  velocity  v  of 
the  plate  is  given  by  solutions  of 


AAv  -  k  v  =  0 


(Ill-i) 


As  in  the  one -dimensional  case,  this  equation  -  and  therefore  the 
impedance  given  below  -  is  only  correct  when  inequality  Eq.  (II-3) 
is  fulfilled. • 


The  impedance  of  point  driven  plates  has  been  calculated  by  sev¬ 
eral  authors--*—*--0 /  in  different  ways.  •  The  result  is 


Zv  =  =  8^/ixn  . 


(III-2) 


In  this  equation  D  is  the  bending  rigidity*  and  m  the  surface 
mass  of  the  plate. 

Equation  (III-2)  also- holds  for  plates  consisting  of  several 
layers,  •  provided  that  they  are  connected  in  such  a  way  that  only 
pure  bending  motion  occurs.  In  this  case  D  is  the  bending 
rigidity  of  the  combination  and  m  its  surface  mass.  For  homo¬ 
geneous  plates  Eq.  (III-2)  can  also  be  expressed  in  terms  of 
Young's  modulus  E  or  in  terms  -of  the  velocity  of  longitudinal 
waves  cT.  Thus  we  get 


*  Because  of  the  lateral  contraction,  the  bending  rigidity  of 
beams  is  always  a  little  smaller  than  the  bending  rigidity 
of  plates.  But.  for  most  practical  cases  they  can  be  set  equal 
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I —  c  ph2 

Zw  =  2.3  h2-\/-^  =2.3  — — rr-  «  2.3  cT  ph2  (lll-2a) 

*  V  .  1-|JL  u 

(p  =  density,  |i  =  Poisson1 s  ratio. ) 

B.  Impedance  of  Infinite  Orthotropic  Plates 

Orthotropic  plates  have  different  bending  stiffnesses  for  vibra¬ 
tions  in  different  directions.  Examples  of  orthotropic  plates 
are:  plates  made  of  non-isotropic  material,  plates  with  grooves 
or  ribs,  corrugated  plates,  or  grillages  consisting  of  crossed 
beams. 

If  the  distance  between  the  grooves,  ribs,  etc.  (whatever  is 
larger)  and  the  plate  thickness  is  much  smaller  than  the  short¬ 
est  bending  wavelength  on  the  plate,  the  bending  vibrations  are 
described  by  the  following  equation  first  obtained  by  Hubert/ : 


D  ^+2D  ----- 

x  a?  ^  3x2Sy2 


a4v  2 

+  D  — tt  —  co  mv  —  0- 

y 


(HI-3) 


Here  D  Is  the  bending  rigidity  in  the  stiff est  direction  (l.e., 
for  example^  in  the  direction  of  the  grooves  or  ribs),  By  is 
the  rigidity  in  the  least  stiff  direction  (i.e.,  perpendicular 
to  the  grooves  or  ribs),  and  m  Is  the  average  surface  mass. 

is  also  a  kind  of  rigidity  which  may  either  be  measured 
(see  the  papers  by  Hoppman^/)  or  calculated. 


For  a  corrugated  plate  (illustrated  in  Fig.  3a)  the  rigidities 
are  given  by1--* 
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D  = 
x 


El 


>  =  £  _ Eh3  .  D*  =  s  fh3 

y  s  12(1_M,2)  xy  g  12(1-hi)  ' 


(HI-*) 


for  a  plate  with  ribs  (illustrated  in  Pig.  3b)  we  have 


Dx  =  EIx  ^  Dy 


In  these  equations  E  is  Young's  modulus,  I  is  the  moment  of  inertia 

A 

along  the  x  axis,  (X  is  Poisson's  ratio,  D'  is  the  torsional,  rigidity 
of  the  plate  without  ribs,  and  B'  is  the  torsional  stiffness  of  one 
rib.  The  dimensions  g,s,h,.u,t,H  are  defined  in  the  figures. 


Finally  for  grills  consisting  of  perpendicular  beams  (illustrated 
in  Fig.'  3c)  we  get 


(III-6) 


Here  B  is  the  bending  stiffness  and  B'  the  torsional  stiffness  of 
the  beams  (see  Fig.  3). 


The  force  impedance  of  orthotropic  plates  for  bending  waves  was 
recently  calculated  by  Heckl.^/  The  result  is 


(III-7) 
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In  this  equation  K(a)  is  the  complete  elliptical  integral  of  the 
first  kind,  for  which  tabulated  values  are  available  (for  example  ' 
see  Ref.  17). 

If  D^r  >  V D  D  the  values  of  a  become  imaginary.  In  this  case 
T  x  y 

the  transformation 

K(ia)  =  — L;  k/  -2= 

■  V^a2  1+a2 

for  a  =  ia  has  to  be  made. 


For  most  practical  cases  D  ~-\/D  D  ,  then  Eq.  (III-7)  can  be 

xy  V  x  y 

approximated  by 


Zp  «  8  (m2  DXD  )1/4  .  (III-8) 

"  «* 

Comparison  with  Eq.  (ill- 2)  shows  that  for  orthotropic  plates 
~\j'. DxCy  is  analogous  to  D  in  isotropic  plates. 


C.  Impedance  of  a  Homogeneous  Plate,  Infinite  in  One  Direction, 
Simply  Supported  in  the  Other  One 

In  this  case  the  impedance  can  be  found  by  expanding  the  velocity 

00 

and  the  force  in  series  of  the  form  ^  f (x)  sin  and  inserting 

these  in  the  bending  wave  equation.  After  some  calculations  the 
impedance  turns  out  to  be  given  by 


1_  =  4i 

ZF  S^Dm 


I 

n=l 


Wn2m2-k2i2"  VnV+k2i2 


I  sin 


nnyo 

l 


(III-9) 
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(£  is  the  distance  between  the  support  edges,  yQ  the  coordinate 
of  the  driving  point).  The  geometry  is  shown  in  Pig.  4. 


As  pointed  out  in  Eq.  (C-20)  of  Appendix  C,  Eq.  (III-9)  can  be 
obtained  also  as  a -limiting  case  of  the  impedance  of  a  beam-plate 
system. 

For  very  small  values  of  £  and  low  frequencies:  i.e.,  for 
k£  <  -Tj  ,  Eq.  (III-9)  can  be  approximated  by 


1__  w  la>£2 
ZF  2tt3D 


sin 


m7yo 

£ 


(HI-9a) 


We  see  that  in  this  case  the  plate  acts  as  a  spring  with  a  stiff¬ 
ness  constant  given  by 


s 


2tt3D 


£  sin 


mryc 

~T~ 


For  very  large  .values  of  £  and  high  frequencies,  i.e.,  for  ki  »  7 r 
the  sums  in  Eq.  (III-9)  can  be  evaluated  to  give  =  8Vl>m  which 
is  the  impedance  of  a  plate  infinite  in  both  directions. 


For  the  intermediate  range,  the  reciprocal  impedance  for  a  plate 
driven  in  the  center  (y  ’  =  -^)  is  plotted  in  Fig.  4  and  Fig.  5* 

The  full  lines  show  the  real  part,  the  dashed  lines  the  imaginary 

part. _ It  can  be  seen  that  for  n  >  9  g  the  average  value  of  , _ v 

Re |  is  very  close  to  unity,  and  the  average  value  of  ImV^— ■  j 
is  nelrly  zero.  ,  Thus  for  high  frequencies  the  impedance  of  suchF 
a  semi-infinite  plate  is  very  close  to  the  impedance  of  an  infinite 
plate  [see  Eq.  (III-2)].  This  result  also  holds  for  plates  finite 
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• 

in  both  directions  (see,  e.g.,  Skudrzyk-^/) .  If  one  is  especially 
interested'  in  the  mechanical  power  transmitted  by  a  point  force 
into  a  finite  isotropic  or  orthotropic  plate,  one  can  show  that 
the-  power  averaged  in  frequency  is  inversely  proportional  to  the 
force  impedance  of  the  infinite  system.  This  fact  can  be  ex¬ 
pressed  by  the  formula 


(III-10) 

The  proof  for  this  equation  can  be  given  in  exactly  the  same  way 
as  it  was  given  in  the  one-dimensional  case. (see  Section  II-A4). 
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IV.  FORCE  IMPEDANCE  OF  BEAM-PLATE  SYSTEMS 


A.  Impedance  of  Systems  Infinite  In  Both  Directions 

This  .problem  (see  Fig.  6a)*  was  solved  by  G.  Lamb^^/  using  the 
notations 

mg  =  mass  per  length  of  the  beam, 
nip  =  surface  mass  of  the  plate, 

B  •=  bending  stiffness  of  the  beam, 

D  =  bending  rigidity  of  the  plate,* 


_!_  _ 

ZF 


ik. 


B 


2TrmBcB 


+00 

/ 


dx 


(kp-x2)  +i(k2+x2) 


^72] 

J 


(iv-i) 


It  should  be  noted  that  nig  and  mp,  and  B  and  D  have  different 
dimensions,  e.g.,  in  the  case  of  rectangular  cross-sections  we 
have 


Egh^b 

=  mp  =  Pp^p  B  =  p2 


Ephp 


12(1 -M.) 


For  dimensions  see  Fig.  6a. 

**  Since  we  assumed  harmonic  motion  of  the  form  e+ia)t  we  have  +i 
where  there  is  -i  in  Lamb's  expressions. 
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For  ^  |  r3  -  kp  ^  »  1 


the  following  approximation  can  be  obtained 


1 

H°b 


1-1-4 


ID 


ir  kgB 


1  + 


3  -  r 
2r3(l-r2) 


cos-^r  + 


A  L 

3  +  r 


o  P  l72 
2r3(l+r2) 


sinh 


4 


(rv-2) 


Physically  Eq.  (IV-2)  shows  that  if  the  beam  is  much  stiff er  than 
the  plate  and  the  frequency  not  too  low,  the  force  impedance  of  a 
beam-plate  system  is  nearly  the  same  as  the  impedance  of  the  beam 
alone  [compare  with  Eq.  (11-4)]. 


On  the  other  hand  it  can  be  shown  that  if  the  stiffness  of  the  beam 
is  not  much  higher  than  that  of  the  plate,  the  impedance  of  the  beam- 
plate  system  is  approximately  the  impedance  of  the  plate  given  by 
Eq.  (III-2). 

B .  Impedance  of  Systems,  Infinite  in  One  Direction  and  Simply 
Supported  in  the  Other  One 

This  problem,  whose  geometry  is  shown  in  Fig.  6b,  is  treated  in 
Appendix  C.  Basically,  the  calculations  show  that  for  beams  much 
stiffer  than  the  plate,  the  beam-plate  system  behaves  like  a  simply 
supported  beam  with  a  loss  factor 


_2  ^ 
kP  “b 


(tv-3) 


and  an  additional  linear  mass 


Report  No.  774 


Bolt  Beranek  and  Newman  Inc. 


AmB  =  4  mP 


(rv-4) 


The  added  damping  and  mass  arise  from  .the  waves  coupled  to  the  plate 
from  the  beam. 
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V.  FORCE  IMPEDANCE  OF  RINGS 

20/ 

According  to  Love — '  the  radial  velocity  v  of  a  ring  with  radius  a, 
linear  mass  m  and  bending  stiffness  B  is  given  by 


B  /d6v 
ma  \d0^ 


+  2 


A  A/\ 
iA  b<p2} 


+  cd2v 


(V-l) 


In  this  equation  pt  is  the  press\ire  in  the  tangential  direction, 
which  in  our  case  is  zero,  p^  is  the  pressure  in  radial  direction, 
which  in  our  case  is  F6(0)  (F  =  driving  force,  6(0)  =  delta  function), 
and  0  is.  the  angle  around  the  ring. 


Equation  (V-l)  is  valid  only  when  the  thickness  of  the  ring  is  much 

smaller  than  one  wavelength  (This  also  means  that  the  thickness 

Ct 

has  to  be  very  small  compared  with  the  radius  . )  and  when  cd  <  — , 

a. 

where  cL  is  the  longitudinal  velocity  in  the  ring  material. 


A  rather  straightforward  method  of  solving  Eq.  (V-l)  under  given 
boundary  conditions  is  to  expand  the  velocity  in  terms  of  the  eigen¬ 
functions  cos  n0.  With  the  use  of  this  method  the  velocity  at  the 
driving  point  may  be  obtained: 


00 


Z  V: 


F  n2 

15555  n^Q  (l+n2)  -"n2(n2-l)aa 


(V-2) 


where  a 


B 

W~2  • 

ma  cjo 


For 


rings  with  rectangular  cross 


section  a  = 


Eh 

~  JT2 

12pa  cd 


(V-3) 
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From  these,  equations  we  get  the  so-called  "modal  impedance" 


F  ioxmra 


JFn  v 


n  n 


l+n2-n2(n2-l)2a 


(V-4) 


Using  the  equation  for  .the  resonance  frequencies  <Dn  of  the  ring 


(see  Love^/ )  • 


2  B  n2(n2-l ) 2 


co„  =  — ^ s 


n 


ma  n  +1 


(y-5) 


we  get 


Zpn  =  icDmira 


1+n 


n 


1  - 


2  , 
.  CD _ 

n 

CD 


(V-6) 


If  the  actual  rather  than  the  modal  impedance  Is  of  Interest,  we 

p 

have  to  carry  out  the  summation  in  Eq.  (V-2) .  For  n  »  1  this 
can  be  done  analytically.  The  result  is 


_!  =  P 

Zp  4ioxna 


ctn  ttB  +  1 


where  p  = 


(1/a) 


1/4 


(V-7) 
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VI.  MOMENT  IMPEDANCE  OF  BEAMS 
A .  Excitation  of  Bending  Waves 

The  moment  Impedance  of  beams*  exclt'ed  Into  flexural  wave  motion 
can  be  found  by  calculations  which  are  nearly  identical  to  those 
used  in  Section  II-A.  This  is  not  surprising  since  the  bending 
wave  equation  for  the  angular  velocity  w  has  the  same  form  as 
Eq.  (II-1). 

The  moment  impedance  of  a  beam  extending  from  x  =  -»  to  x  +«  driven 
at  x  =  0  is 


z„  _  2E§A  (1+i)  .  ZF 


M 


0) 


k 


2 


(VI-1 ) 


[Zp  =  force  impedance  given  by  Eq.  (II-4) ] 


For  a  beam  extending  from  x  =  C  to  x  =  +°°  driven  at  x  =  0,  we  get 


(under  the  boundary  condition 


d2w 


dx 


=  0) 
x=0 


z  _  h. 

M_  k2 


(VI-2) 


[Zp  =  force  impedance  given  by  Eq.  (II-5)]. 


*  For  the  geometries  considered  here  see  Fig.  1  and  replace  the 
force  by  a  moment,  e.g.,  by  two  opposite  forces  a  short  distance 
apart . 
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It  is _now  very  easy  to  find  the  moment  impedance  o&  beams  with 
reflecting  devices;  all  we  have  to  do  is  to  use  the  equation  in 

p 

Section  II,  to  replace  me  by  mc/k  ,  and  to  insert  the  reflection 
coefficients  rw  for  angular  velocities  instead  of  the  reflection 
coefficient  r  for  flexural  velocities.  The  latter  substitution  is 
no  problem  at  all,  because  it  turns  out  that  in  general  r  =  -r. 
This  means  that  the  list  in  Appendix  B  is  useful  for  moment  im¬ 
pedances,  too*. 

Thus  we  can  summarize:  The  moment  impedance  of  a  beam  can  be  found 

2 

from  the  force  impedance  of  the  same  beam  by  replacing  me  by  mc/k 
and  r  by  -r. 


B .  Excitation  of  Torsional  Waves 

If  a  beam  is  excited  by  a  moment  in  such  a  way  that  torsional  waves 
are  excited,  the  moment  impedance  can  be  found  by  solving  the  tor¬ 
sional  wave  equation 

w  =  0  (VI-3) 

dx  CT 

for  the  proper  boundary  conditions**.  In  Eq.  (VI-3) 


Cip 


(VI -4) 


*  The  reflection  coefficient  Rw  for  the  near  field  could  also  be 
found  in  a  very  simple  way  because  of  the  general  relation 
Rw  =  iR. 

**  Strictly  speaking  Eq.  (VI-3)  is  true  only  for  cylindrical  beams, 
but  it  is  also  a  very  good  approximation  for  rectangular  beams 
as  long  as  the  ratio  of  the  lengths  of  the  two  sides  is  not  too 
large  compared  to  unity . 
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is  the  propagation  velocity  for  torsional  waves  and  G  =  g (1+pT)  is 
the  shear  modulus . 

Since  Eqs .  (VI-3)  and  (VI-4)  are  very  'similar  to  Eqs.  (11-23)  and 
(11-24),  the  moment  impedances  for  torsional  waves  are  also  very 
similar  to  the  force  impedances  for  longitudinal  waves. 


For  a  beam  extending  from  x  =  +»  to  x  =  -»,  driven  somewhere  in 
between  we  get 


(VI -5) 


(S  =  cross-section  of  the 


beam,  k  =  radius  of  gyration) . 


For  a  beam  extending  from  x  =  0  to  x  =  »,  driven  at  x  =  0  we  get 


ZM  =  SK2y~Gp" 


(VI-6) 


For  beams  with  reflecting  devices  we  can  use  Eq.  (lI-27)-(lI-29), 

p. 

if  we  replace  S  by  Sk  and  E  by  G. 


i 
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VII.  MOMENT  IMPEDANCE  OP  PLATES 

In  the  case  of  a  force  act-ing  on  a  plate  it  is  possible  to  give  an 
impedance  which  is  independent  of  the  shape  and  size  of  the  area 
over  which  the  force  is  acting,  provided  that  the  area  is  suffi- 
ciently  small.  Unfortunately  the  same  does  not  hold  in  the  case 
of  moment  impedance. 


If  the  moment  is  applied  in  form  of  a  dipole  of  forces  with  mutual 
distance. d,  the  moment  impedance  found  by  using  Eq.  (III-1)  is, 
according  to  Cremer^/, 


Z 

^M  “  CD  1 


i  +  j:  (l-ln2kd) 


(kd  «  1) 


(VII-1) 


On  the  other  hand  for  a  moment  applied  over  a  disc  of  radius  a 
Dyei^/  found  using  Eq.  (ill-l) 


= 


4d 


JM  "  O)  1  i  .  . 

¥  -  ?  ln  ka 


(ka  «  1) 


(VII-2) 


We  see  that  the  moment  impedance  is  not  independent  of  the  shape 
and  area  over  which  it  is  applied;  furthermore  for  a— ^  0  the 
absolute  value  of  the  impedance  would  go  to  zero. 


These . difficulties  still  exist  if  the  more  accurate  bending  wave 
equations  by  Mindlin^/  are  used;  these  equations  also  take  into 
account  the  local  shear  deformations  which  always  exist  near  the 
source.  In  this  case  Dyer^/  found  for  the  impedance  applied  over 
a  disc  of  radius  a 
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ZM  =  ~ 


2D( 1+L) 


cc 


1 

% 


I  In  ka  + 
t r 


1  2L_  /_hf ' 1 
7T  l-|i  l  ?raj  J 


(VII-3) 


In  this  equation  h  is  -the  thickness  -of  the  plate .  The  quantities 
L  and  the  third  term  in  the  denominator  can  be  found  from  Fig.  J. 

Although  Eqs .  (VII-1 )-(VII-3)  give  very  different  values,  they 
have  one  common  feature --the  real  parts  in  the  denominator  are 
equal.  This  fact  is  of  physical  significance  since  the  power  trans¬ 
mitted  by  a  localized  moment  is 


P  =  ||M2|Re|J 


M 


(VII-4) 


In  the  first  two  cases  this  gives 


P  = 


M2  co 
3  2D 


M 


8  \j  Dm 


(VII-5) 


in  the  third  case  we  get 


p  _ 

"  l6D(l+L) 


(VII-5a) 


which  is  not  very  different  from  Eq.  (VII-5). 

In  summary:  If  only  the  mechanical  power  transmitted  by  a  localized 
moment  is  of  interest,  the  area  and  shape  over  which  the  moment  is 
applied  is  not  of  interest.  In  this  case  the  power  is  approximately 
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given  by  Eq.  (VII-5) .  If*  however,  the  impedance  itself  is  important 
(e.g.,  when  it  has  to  be  added  to  a  source  impedance),  it  is 
necessary  to  use  Eq.  (VII-3) • 
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VIII.  MOMENT  IMPEDANCE  OF  RINGS 


The  moment  Impedance  of.  rings  can  be  found  in  a  way  that  is  very 
similar  to  that  used  .in’  Section  V.  If  a  dipole  of  forces  with 
moment  M  is  acting  in  the  radial  direction,  the  angular  velocity 
at  the  excitation  point  is 


oo  h 

_ M _ 1  yi  _ n_ _ 

ltuOTa  ( 2ira ) 2  ^  (l+n2)-n2(n2-l)2a 


(VIII-1) 


with 


a  = 


B 

jj— ? 

ma  co 


From  Eq.  (VIII-1 )  we  find  for  the  modal  moment  Impedance 


JMn 


=  iourmra 


l+n‘ 


n 


0) 


1- 


n 


O) 


(2jra 
n  , 


(VIII -2) 


(con  see  Eq.  V-5) . 

p 

The  asymptotic  value  of  l/Z^  for  n  »  1  turns  out  to  be 


^M  l67T^icoma^ 


where  p  =  [l/a'j 


1/4 


(ctn7rp  +1) 


(VIII-3) 
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APPENDIX  A 

LIST  OP  REFLECTION  COEFFICIENTS  FOR  BENDING  WAVES  ON  BEAMS* 

•1 .  Change  In  Cross-Section  or  Material 

If  m^,  B^,  denote  the  mass,  bending  stiffness,  and  wave  number 
before  the  discontinuity,  and  mg,  Bg,  kg  the  same  parameters  behind 
it,  the  reflection  coefficient  r  for  the  wave  field  is  given  by  (see 
Cremer^/ ) 


r  =  2a(i-P2)  -  iP(l-q)2 
p(l+q)2  +  2cx(l+p2) 


(Al) 


For  the  near-field  reflection  coefficient  Cremer  obtains 


r'=  P(l-q2)-  iS(l~q2) 

P (l+a) 2  +  2a(l+p2) 


(A2) 


where 


(flip®?  ^2 

a  and  p 


2.  Impeding  Masses 

If  mQ,  0,  K  denote  the  mass,  rotatory  inertia  and  the. radius  of 
gyration  of  the  added  mass,  and  m,  k  the  surface  mass  and  bending 
wave  number  of  the  beam,  the  reflection  coefficients  are  given  by 
(see  Cremer^ ) 


*  For  sketches  of  the  geometries  see  Fig.  8. 
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and 


r  = 


2  n  2  I 
-M.+eV+-|e  M- 


( LL+  e  V )  -  i  (  4+(.l-  e  2(i  3--i  £  V1 ) 


(A3) 


R  = 


-  e  2|X3- V+l  ( 2M-4 ) 
(|i+€2p,-3)  -l(4+p.-  e2(jL3-4e2ti4) 


(A4) 


m 

Here  the  abbreviations  £  =?  —  K  and  [X  =  —  k  are  Introduced. 

m  m 

o 

For  frequencies  more  than  one  octave  above  the  "Sperrfrequenz"  which 
Is  given  by  kH  =  1  the  above  equations  can  be  approximately  written  as 


i-M2  *  .(^)2  ;  R  *  -(i-i)  •  ;  (A5) 

3 .  Corners  and  Other  Junctions' 

a)  Vibrations  in  the  Plane  Formed  by  the  Beams . 

This  case  was  investigated  by  Cremer^  taking  into • consideration  the 
excitation  of  longitudinal  waves.  The  corresponding  formulas  can  be 
found  in  Cremer's  paper.  For  many  practical  cases,  however,  it  is  not 
necessary  to  use  these  rather  complicated  formulas,  since  the  errors 
are  not  too  great  in  the  region  where  the  simple  bending  wave  theory 
is  valid  [see  Eq.  ( II— 1 ) ] .  Therefore,  for  most  cases  the  approximate 
formulas 


r 


1+iq 

1+q 


R  =  - 


(A6) 


are  sufficient. 


-35- 


Report  No.  774 


Bolt  Beranek  and  Newman  Inc. 


Here  q  is  given  by  the.  bending  stiffness. B  and  the  wave  number  k 
of  the  primary  beam  and  the  same -quantities  Bv  and  kv  of  the  • 
adjoined  beams 


«  -  IE  E 

v=l 


B  k 


v  v 


(A7) 


b .  Vibrations  Perpendicular  to  the  Plane  Formed?  by  the  Beams 

If  this  case  is  treated  accurately ,  . the  excitation  of  torsional 
waves  must  be  considered.*  This  leads  to  the  rather  complicated 
expressions  given  in 'Ref.  11.  If  the  .excitation  of  torsional  waves 
is  neglected  the  reflection  coefficients  are- 


3 


(A8  )* 


and 


R  =  q  •  (A9) 

The  quantity  q  Is  given  by  Eq.  (A7)  but  since  the  bending  stiffness 
is  in  general  different  along  different  axes  the  numerical  value  of 
q  might  be  very  different . 


Equations  (A8)  and  (A9)  are  not  very  good  approximations,  and  are 
especially  poor  at  low  frequencies  .  They  can  serve  'only  as  very 
rough  estimates . 


Terminations 

For  clamped  ends: 

V 

r  =  -i; 

R  =  -(1-i)  . 

lA10) 

For  simply  supported 

ends:  r  =  -1 ; 

R  =  0  . 

(All) 

For  free  ends: 

r  =  -i; 

R  =  1-i  . 

(A12) 
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If  it  is  possible  to  describe  the  termination  by  a  force  impedance 
P  M 

—  =  Zp  and  a  moment  impedance  —  =  Z^,  the  reflection  coefficients 
are  given  by 


and 


(A13) 


R  -  "21  ( 1-TV )  .  /  n-l  J,  \ 

R  “  •(1-1)(1+tv) -2(r-iv)  J  •  (Al4' 

where  v  -  -  |g  Z„  and  r  =  ^  Zp  . 

The  reflection  coefficients  for  free,  clamped  and  supported  beams 
can  easily  be  verified  from  these  equations.  Also,  the  case  of  a 
beam  attached  to  a  plate  or  to  a  large  mass"ls"  included  in  Eq.  (A13) 
and  (Al'4) .  It  should  be  mentioned  that  |r|  is  always  unity  when  v 
and  t  are  real'i.e.,  when  Zp  and  are  purely  imaginary. 
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APPENDIX  B 

IMPEDANCE  OF  BEAMS  WITH  REFLECTING  DEVICES 

A  model  of  a  beam  with  two  reflecting  devices  is  shown  in  Fig.  le. 
To  derive  the  -velocity  at  the  driving  point  we  assume  that  L  and  l 
are  at  least  larger  than  half  a  bending  wavelength;  this  implies 
that  the  near  fields  excited  at  the  discontinuities  can  be  ne¬ 
glected  . 

The  velocity  at  the  driving  point  x  = -0  can  be  calculated  by  noting 
that  the  primary  waves  excited  at  x  =  0  propagate  to  x  =  I  or  . 
x  =  -L.  At  these  points  they  are  partially  reflected  and  propagate 
.in  the  opposite  direction  until 'they  encounter  the  next  disconti¬ 
nuity  where  they  are  reflected  once  more,  etc.  The  addition  of  all 
these -partial  waves  gives  the  actual  velocity. 

The  velocity  of  the  primary  waves  is  given  by  Cremer^  as 


and 


-ie  kx)  for  x  ^  0 


v 


F  t  +ikx  ,  kx\ 

=  TfHST  (e  "ie  } 


for.  x  ^  0 


(Bl) 


If  we  first  consider  the  v+  wave  only,  we  find  (neglecting  near  fields) 
that  after  the  first  reflection  the  wave  has  ,the  form 


v 


R1 


F 

4cm  rj 


s-ik(2  i-x) 


(B2) 


(r^  =  reflection  coefficient  at  x  =  l). 
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If  the  wave  described  by  Eq.  (B2)  hits  the  discontinuity  at  x  =  -L, 
another  reflected  wave  is  generated,  namely 


v  =  _2i  rflr  e"ik(2i+2L+x) 

VR2  Ifcm  ri.L  e 


*(B3) 


The  next  reflections  at  x  =  &  and  x  =  -L  give 


and  • 


(B4) 


v  ^  r  r  e-ik(2^+2L+x)  *  -2ik(L+A) 

VR4  7cm ‘r£rL  e  rL Vl  6 


and-  so  on.  I-t  can  be  seen  immediately  that  the  summation  of  all  of 
the  reflections  forms  a  geometrical  series  .in-  r^r^  e_2ik(L+i!).  The 
summation  therefore  can  be  carried  out  easily,  and  if  we  do  this  for 


v+  and  v_  we  obtain 


F  l-i+rL  e  21kL+rg  e“2iki+(l+i ): 
V°  =  ^  i-r^L  e"2lk^+: 


+(l+i)rgrL  e~2ik(i+L) 


(v  =  velocity  at  x  =  0). 
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APPENDIX  C 

FORCE  IMPEDANCE  OF  A  FINITE  PLATE  BEAM  SYSTEM 

The  input  impedance  of  an  infinite  bar  which  is  attached  over  its 
entire  length  to  an  infinite  plate  was  derived  by  Lamb^/,  An 
extension  of  this  problem  will  be  treated.,  namely,  the  force  im¬ 
pedance  of  a  finite  beam  attached  to  a  plate.  The  geometry  of 
the  system  can  be  seen  in  Fig.  6b.  To  simplify  the  calculations 
it  is  assumed  that  the  plate  and  the  beam  are  simply  supported 
at  the  lines  y  =  0  and  y  =  i,  i.e.,  the  velocities  and  the  mo¬ 
ments  at  these  lines  are'. .zero. 

We  assume  the  plate  and  the  beam  to  be  thin,  so  that  the  vibra¬ 
tions  of  the  plate  are  given  by 


AAvp-  “  kp  vp  =  0 

and  the  vibrations  of  the  be  Am  by 

" 


,4 

d  v. 


B 


dy 


*5 


VB  “ 


ioo 

B  PB 


(Cl) 


(C2) 


In  these  and  the  following  equations  the  subscript  B  corresponds  to 

the  bending  wave  numbers,  B  the  bending  stiffness  of  the  beam,  and 
D  the  bending  rigidity  of  the  plate.*  Pg  is  the  pressure  acting  on 


the  beam  and  P  to  the  plate  I  A  is  the  Laplacian  operator,  kg 


*  See  note  on  p.  22. 
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the  beam  and  consists  of  two  parts,  the  exciting  pressure  pBE  given 
by  a  point  force  P  applied  at  a  point  y  .  and  a  pressure  p™  which 

O  Jda 

includes  the  backward  reaction  from  the  plate. 


Because  of  the  simple  boundary  conditions,  the  velocity  of  the 
beam  can  be  expressed  as  . 


00 


n=l 


(C3) 


To 'find  the  unknown  coefficients  vBn;  .we  must  know  the  pressures 
PB  and  PBpj.  To  this  end  we  expand  the  velocity  of  the  plate  in  a 
series- similar  to  Eq.  (C3).  (This  can  be  done  because  of  the 
similarity  of  the  boundary  conditions.)  We  write 


VP  .=  Z  fn(x)  sin  .  •  (C4) 

n=l  .  • 

The  functions  fn(x)  describe  the  propagation  o'f  bending  waves  into 
the  plate.  They  can  be  found  by  inserting  Eq;  (C4)  into  Eq.  (Cl) 
giving  •  - 


,4  2  2 

-TF  f  (*)  -2 
•  dx  &  dx 


5  fn<x>+ 


f  (x) 
n v  ' 


(C5) 


This  equation  can  easily  be  solved  by  introducing  the  exponential 
function.  Thus  we  get 
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f  (x)  -  v  e1f>x+v '  eax+v  e’1Px+v'  e"0^ 
*nw  "  vPn-e  +vPn-e  +vPn+e  +vPn+e 


where 


P  =  \/k2  -^and  cAhl  +  *2*2 


^  i2 


(c6) 


(C7) 


For  our  model  we  must  consider  that  no  waves  can  come  from  infinity 
Therefore,  we  have  two  different  solutions  for  x  <  0,  and  x  >  0. 


P+ 


-E 


vPn+e'ipX+vPn+e’°X  )  sln  T  »  for  x  >  0 


n=l 


and 


(C8) 


=  Z  f  Pn-elpX«Pn-eC,X  )  sln  y  for  x  <  0 


n=l 


Because  of  the  symmetry  of  the  problem  we  can  put  vpn+  =  vpn_ 

and  v '  .  =  v '  .  Furthermore,  it  is  sufficient  to  consider  one 

rn+  Fn-*  _  _ 

of  the  equations  (C8)  and  to  double  the  reacting  pressure  obtained 
in  this  way . 


The  coefficients  vpn+  and  v^n+ 
by  using  the  following  boundary 


can  be  expressed,  in  terms  of  vBn 
conditions . 


VB  "  vP+j 


or  v 


Bn 


=  vpn+  for  x  =  0 


(C9) 
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and 

(CIO) 

The  first  of  these  equations  expresses  the  fact  that  the  velocity 
of  the  beam  and  the  plate  must  be  equal  at  the  line  x  =  0.  The 
second  equation  means  that  the  angular  velocity,  at  x  =  0  must  be 
zero,  because  there  are  no  bending  moments  in  the  x  direction. 

Equation  (CIO),  together  with  Eq.  (C8)  and  Eq„  (C9),  gives 

(Cll) 


(C12) 


v. 


Pn+ 


Bn  (c2+ipo) 


2k: 


and 


v. 


•v 1 

Pn+ 


Bn  /  2 


2k: 


(p  -ipa) 


dv 


P+ 


*55r 


0  for  x  =  0 


Now  we  are  able  to  express  the  pressure  pBR  in  terms  of  vRn.  The 
pressure  is  given  by 

(C13) 

i  =  0 

Equations  (C8),  (Cll),  (C12),  and  (C13)  yield,  after  some  algebra, 

00 

pBR='ll  Z  vBn(p2a  ipa2)  sin  ^  y  .  (Cl4) 

n=l 


_  2D 
PBR  "  ico 


dAv 


P+ 


dx 
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Now  we  can  expand  the  exciting  pressure  pBE  in  the  same  series  by 
writing 


Pbe 


-  E 


BEn 


n=l 


(C15) 


where  PBgn  can  easily  be  computed  in  the  special  case. of  a  point 
force  '  • 


%E  =  PS(y-s0) 


(Cl6) 


which  has  the  Fourier  coefficients 


mry 

PBEn  =  I  F  Sin  ~~i 


.ciT) 


By  inserting  Eqs .  (C13),  (Cl4),  (C15),  and  (C17)  in'Eq.  (C2)  we 
get 


mr 1 
1  l  I 


2  ico  ^0 

v„  =  4  F  sin  - - 

Bn  i  B  i 


+  2  js(op2-ip?2)vBn 


(018) 


This  expression  gives  us  the  quantities  vBn-  We  get  finally 


,  _  2icbF  V 

B  ~Br~  ^ 


<  n7ryo  ,  mr 
sin  — j —  sin  -j-  y 


-if2?)  ->4-a  ife-  -4 )L 


(C19) 
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This  equation  is  very  similar  to  that  derived  by  Lamb^/ ;  the  main 
difference  is  that  a  sum  instead  of  an  integral  is  obtained  for 
the  finite  case. 

To  obtain  some  further  information  from  Eq.  (C19)  let  us  investi¬ 
gate  two  special  cases:,  first,  the  case'of  a  very  light  and 
flexible  beam  on  a  stiff  plate;  second,  the  case  of  a  very  stiff 
beam  attached  on  a  rather  flexible  plate; 


Case  A 

In  this  case,  especially  when  there  is  no  beam,  i.e.,  B  =  0,  we 
get,  for  the  velocity  at  the  driving  point 


'BO 


~  + 


tuF 
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E  n=l 


sin 


2  n^o 


_  2  » 2  2  2 
Ly^p  a  — n  tt 


4 


'kp  i2  +  n2ir 2 


.  (C20) 


Equation  (C20)  is  identical  with  Eq.  (III-9),  as  expected. 


Case  B 

In  this  case,  which  occurs  very  often  in  practice  (e.g.,  when  a 

thin  plate  Is  stiffened  by  a  stringer),  it  is  reasonable  to  compare 

the  vibrations  of  the  beam-plate  system  with  the  vibrations  of  a 

damped  beam. 

•  ■ 

In  the  case  of  a  point -driven  damped  beam  we  insert  the  complex 
value”" 

=  *4  (1-m) 
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In  Eq.  (C2)..  The  vibrations  can  be  calculated  as  shown  above. 
The  only  difference  is  that  there  is  no  reacting  pressure  PBR- 
Thus  we  get  for  the  velocity  of  a  point -driven  damped  beam 


00 

!icnP  V 
B T  L 


sin 

(rnr 


mry 


T 


,  4  .  ,  4 


-sin 


(C22) 


To  compare  this  equation  with  Eq.  (C19)  we  take  into  consideration 
that  for  D  «  Bkp  the  value  of  the  sum  in  Eq.  (C19)  is  mainly 
determined  by  the- behavior  .of  the  denominator  in  the  vicinity  of 
rnr/i  =  kg.  Therefore,  we  make  no  great  error  if  we  put  n ir/£  =  kg 
in  the  last  term  of  the  denominator  of’Eq.  (C19) •  So  we  get 
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If  we  compare  this  result  with  Eq.  (C22),  we  see  that  the  coupling 
of  a  flexible  plate  to  a  stiff  beam  has  approximately  the  same 
effect  as  damping  the  beam  with  a  loss  factor 


(C24) 


and  changing  the  bending  wave  number  because  of  an  "additional  mass 
of  magnitude 
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(C25) 


These  two  equations  can  he  simplified  a  little  more  if  we  set 

4  4  » 

.  kp  »  kg,  which  is  at  least  correct  if  beam  and  plate  are  made 

out  of  the  same  material .  Hence 
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j  L 

X  =  0 

(a)  eq.  (H,4)  (b)  eq.  (II,  5)- 


u  i 

X  =  -i  X  =  0  .  '  •  x=0  • 

.•  . 

(c)  eq.  (11,6)  (d)  eq:  (H  J) 


r 


X  =  i 


F 


LJ  |  U 

X  = - L  X=i 

(e)  eq.  (n,  10  AND  II) 


NOTE: 

TO  FIND  THE  MOMENT  IMPEDANCE 
USE  THE  CORRESPONDING 
EQUATION  FOR  FORCE  IMPEDANCE 
AND  REPLACE  cm  BY  cm/k2 
AND  r  BY  -r. 


FIG.  1 


LIST  OF  BEAM  IMPEDANCE  (BENDING) 
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\ a )  eq.  (n,  25) 


(b)  eq.  (E,  26) 


(c)  eq.  (11,27) 
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BEAM  IMPEDANCES  (LONGITUDINAL) 
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(a)  CORRUGATED  PLATE 


i  ~i  n  n  n  n  n  i  r 


(b)  PLATE  WITH  RIB-6 


BEAMS  OF  STIFFNESS  B2  AND  D 


BEAMS  OF-STIFFNESS 
B,  AND  D,‘ 


(C)  GRILLS 


FIG.  3  DIMENSIONS-  OF  ORTHOTROPIC  PLATES 

(impedances  eq.  HI, 4-111,7) 


FIG.4  F 


FIG.5  FORCE  ADMITTANCE  OF  A  PLATE  INFINITE  IN  ONE  DIRECTION,  SIMPLE 
SUPPORTED  IN  THE  OTHER  ONE.  (DRIVING  POINT  IN  THE  MIDDLE) 


CROSS-SECTION,  SEE  UPPER  SKETCH 

(b)  FINITE  SYSTEM 


FIG. 6  SKETCH  OF  BEAM  PLATE  SYSTEM 
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FIG.  7  QUANTITIES  FOR  THE  CALCULATION  OF  THE 
THIN  PLATE  MOMENT  IMPEDANCE 
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